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ORIGINAL ARTICLE 

Investigation of manic and euthymic episodes identifies 
state- and trait-specific gene expression and STABl as a new 
candidate gene for bipolar disorder 

SH Witt\ D Juraeva^ C Sticht^ J Strohmaier\ S Meier^ J Treutlein\ H Dukal\ J Frank\ M Lang\ M Deuschle^ TG Schulze^ 
F Degenhardt^'^ M Mattheisen^'^^ B Brors^ S Cichon", MIVl Nothen^'^ CC Witt^^ and M Rietscher 

Bipolar disorder (BD) is a highly heritable psychiatric disease characterized by recurrent episodes of mania and depression. To 
identify new BD genes and pathways, the present study employed a three-step approach. First, gene-expression profiles of BD 
patients were assessed during both a manic and an euthymic phase. These profiles were compared intra-individually and with the 
gene-expression profiles of controls. Second, those differentially expressed genes that were considered potential trait markers of 
BD were validated using data from the Psychiatric Genomics Consortiums' genome-wide association study (GWAS) of BD. Third, the 
implicated molecular mechanisms were investigated using pathway analytical methods. In the present patients, this novel 
approach identified: (i) sets of differentially expressed genes specific to mania and euthymia; and (ii) a set of differentially expressed 
genes that were common to both mood states. In the GWAS data integration analysis, one gene {STABl) remained significant 
(P= 1.9x 10""^) after adjustment for multiple testing. STABl is located in close proximity to PBMRl and the NEK4-ITIH1-ITIH3-ITIH4 
region, which are the top findings from GWAS meta-analyses of mood disorder, and a combined BD and schizophrenia data set. 
Pathway analyses in the mania versus control comparison revealed three distinct clusters of pathways tagging molecular 
mechanisms implicated in BD, for example, energy metabolism, inflammation and the ubiquitin proteasome system. The present 
findings suggest that STABl is a new and highly promising candidate gene in this region. The combining of gene expression and 
GWAS data may provide valuable insights into the biological mechanisms of BD. 

Translational Psychiatry (2014) 4, e426; doi:1 0.1 038/tp.201 4.71; published online 19 August 2014 



INTRODUCTION 

Bipolar disorder (BD) is a complex and highly heritable psychiatric 
disorder with a lifetime prevalence of 0.5-1.5%, and ranks among 
the top 10 leading causes of years lost because of disability 
worldwide.^ BD is characterized by the occurrence of recurrent 
episodes of mania and depression. Mania is a state of elevated 
mood, which is accompanied by a variety of distinctive symptoms 
such as increased activity, reduced need for sleep, exaggerated 
estimation of ability, rapid thoughts and speech, and irritability.^ In 
contrast, the depressive phases of BD are characterized by 
sadness, anhedonia and fatigue. In between these episodes, the 
mood state is described as euthymic, that is, diagnostic criteria for 
mania and depression are not met. However, research suggests 
that subclinical mood symptoms may persist during euthymic 
periods.^ 

At the time of writing, the etiology and underlying biological 
mechanisms of BD remain largely unknown, and this lack of 
knowledge limits the development of more effective treatments, 
objective diagnostic tools and biological markers of disease 



susceptibility. However, genome-wide association study (GWAS) 
analyses have identified novel genetic risk variants for BD that had 
not been previously implicated in the disorder, and thus provided 
promising new leads into the biology of BD.^'^ Polygenic score 
analyses of these data have shown that many more contributing 
variants remain undetected, as most do not pass the threshold of 
genome-wide significance because of multiple testing.^ Further- 
more, gene-expression studies and state-of-the-art analytical 
strategies, such as pathway- and gene-network analyses, have 
generated new insights into the etiology of BD through the 
identification of novel biological mechanisms.^ To obtain addi- 
tional information, various authors have advocated a convergent 
approach, involving the combining of data from different 
sources.^"^° 

The aim of the present feasibility study was to apply a novel 
three-step convergent approach to identify new candidate genes 
and pathways for BD. First, the gene-expression profiles of a 
sample of BD patients were analyzed during both a manic phase 
(BD manic) and an euthymic phase (BD euthymic). These profiles 
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were first compared on an intrapatient basis, and then with the 
gene-expression profiles of controls. Second, those differentially 
expressed genes from step one that were considered potential 
trait markers of BD were integrated with existing GWAS data. 
Third, the implicated molecular mechanisms were investigated 
using pathway analytical methods. 

MATERIALS AND METHODS 

Subjects 

All patients and controls were of German descent. To enhance sample 
comparability, all blood samples were drawn in the morning during 
routine clinical phlebotomy. The study was approved by the respective 
local ethics committees. Written informed consent was obtained from all 
subjects following a detailed explanation of the study. For patients 
recruited in a manic phase who lacked the capacity to provide informed 
consent, preliminary consent was obtained before blood sampling by the 
treating psychiatrist. The patient was then interviewed after recovery by a 
member of the study team, and was provided with a detailed explanation 
of the study and asked to sign the full informed consent documentation. 
Blood samples were only processed after full informed consent had been 
obtained. In cases where this consent was refused, the respective blood 
samples were destroyed. All subjects included in the present study had 
participated in a previous GWAS of BD (for details see Cichon et al.^). 

BD patients. Blood from both a manic and an euthymic phase was 
collected for 95 BD in-patients. The diagnosis was assigned using the 
Structured Clinical Interview for the DSM-IV axis I disorders (SCID-I)," and 
the Operational Criteria Checklist for Psychotic Illness (OPCRIT v3.32),^^ and 
through the evaluation of medical records and family history. Mania was 
evaluated using the Young Mania Rating Scale (YMRS), a diagnostic 
questionnaire used to evaluate manic symptoms in patients.^^ Details of 
the medications prescribed at each blood sampling time point are 
indicated in Table 1. No statistically significant differences in terms of the 
type of medication prescribed were observed between mood states. 

Selection of BD patients for the gene-expression analyses. An extreme 
group design using male patients was applied. This approach was selected 
as the detection of differential gene expression is optimized by 
comparisons between two groups that show extreme differences in terms 
of the phenotype of interest, but are otherwise homogeneous. Thus, 
patients with the highest differences in YMRS score between manic and 
euthymic states (mean difference: 20.09 ±7.3) were selected (n = 1 1; mean 



age = 48.27 ± 1 0.29 years. Table 1 ). Patients had a mean YMRS score of 25.2 
(±10.2) when manic, and 4.3 (±4.7) when euthymic (Table 1). 

Controls for the gene-expression analyses. Male, age-matched healthy 
controls (mean age = 47.6 ± 1 1.98 years) were recruited for the gene- 
expression analyses (Table 1). 

Blood sampling, RNA preparation and microarray procedure 
Blood samples were collected in RNA-stabilizing PAXgene tubes (Pre- 
Analytix, BD, Heidelberg, Germany), and stored at -80°C until analysis. 
RNA was isolated with the PAXgene Blood RNA Kit 50 (PreAnalytix, Qiagen, 
Hilden, Germany) according to standard protocols. Total RNA yield was 
determined using the Quant-iT RiboGreen RNA Reagent and Kit (Life 
Technologies, Darmstadt, Germany), and a Wallac Victor 2 1420 Multilabel 
Counter (Perkin Elmer, Rodgau, Germany). Total RNA purity was evaluated 
via OD measurements (260 nm/280 nm) in a NanoDrop (peqLab, Eriangen, 
Germany). RNA integrity was determined by RNA integrity number 
measurement using RNA 6000 Nano Assay RNA chips run in an Agilent 
2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The 
inclusion criteria were as follows: a ratio of 1.9-2.2 (OD 260/280); RNA 
integrity number > 8.0, and the absence of a peak of genomic DNA 
contamination in electropherograms. Samples were transcribed to cDNA 
and hybridized to an Affymetrix Gene Chip Human Exon 1.0 ST Array 
(Affymetrix, Santa Clara, CA, USA) using the Whole Transcript Sense Target 
Labeling Assay protocol and lOOng of total RNA from each sample. 

Microarray data analyses 

Gene-expression profiling was performed using Affymetrix Human Exon 
1.0 ST arrays. Biotinylated antisense cRNA was prepared according to the 
Affymetrix standard labeling protocol. One control sample performed 
poorly and was excluded from the analyses. Hybridization on the chip was 
performed in a GeneChip Hybridization oven 640. Chips were dyed in a 
GeneChip Fluidics Station 450 and scanned with a GeneChip Scanner 3000. 
All equipments were manufactured by Affymetrix (Affymetrix, UK). Data 
importing was performed using the Affymetrix Expression Console 
Software package and was restricted to core probe sets that correspond 
to genes. The core set on the raw fluorescence intensity values was 
normalized by applying robust multi-array average. Differential gene 
expression was analyzed using analysis of variance and the commercial 
software package SAS JMP7 Genomics, V4 (SAS Institute, Cary, NC, USA). A 
false-positive rate of a = 0.05 with false discovery rate correction was used 
as the level of significance (threshold: P < 2.59x 10~^). All raw and 
normalized data have been deposited in the Gene-expression Omnibus 
database (http://www.ncbi.nlm.nih.gov/geo/; accession no: GSE46416). 

The following gene-expression profile comparisons were performed: (1) 
BD euthymic versus controls; (2) BD manic versus controls; and (3) BD 
euthymic versus BD manic. 

Validation using qRT-PCR 

cDNA was generated per Life Technologies High-Capacity cDNA Reverse 
Transcription Kit. Triple qRT-PCR reactions per sample were performed 
using Taqman gene expression probes (Life Technologies) for stabilini/ 
STABl (Hs01109068_m1/Exon 25-26) by the 7900HT QRT-PCR system 
(Life Technologies). Fast program: 50 °C, 2:00 min; 95 °C, 0:20 min, 
and 40 repeats of 95 °C, 0:01 min and 60 °C, 0:20 min. P-Actin/ACTB 
(Hs9999903_m1) normalized the relative expression of the STABl gene. 
Fold change differences were calculated by the A Q method. 

Integration of GWAS data 

To explore whether the genes that were differentially expressed in the 
present BD patients and controls were significantly enriched for 
association in patients and controls from a GWAS of BD, gene-based 
association tests were performed on the Psychiatric Genomics Consortium 
(PGC)-BD primary data set.^ The gene-based tests were calculated using 
VEGAS, which accounts for gene size, single-nucleotide polymorphism 
(SNP) density and linkage disequilibrium between SNPs.^^ To capture 
regulatory regions and SNPs in linkage disequilibrium, the default setting 
in VEGAS was used. 



Table 1. Demographics of patients and controls, and clinical 
characteristics and prescribed medication, at the time of blood 
sampling for patients 




Euthymia 


Mania 


Controls 


Group 










differences 


Demographics 










N 




11 


10 




Age, years 


48.3 


±12.0 


47.6 ±10.3 


n.s. 


Clinical characteristics 










YMRS Score 


4.3 ±4.7 


25.2 ±10.2 




0.00 (0.01) 


Medication n (%) 










Antipsychotics all 


11 (100) 


11 (100) 




n.s. 


Typical 


2 (18) 


5 (45) 




n.s. (0.021) 


antipsychotics 










Atypical 


10 (91) 


10 (91) 




n.s. 


antipsychotics 










Mood stabilizers 


4(36) 


4(36) 




n.s. 


Anticonvulsants 


9(82) 


10 (91) 




n.s. 


Antidepressants 


1 (9) 


0(0) 




n.s. (0.038) 


Benzodiazepines 


5(45) 


5 (45) 




n.s. 


Abbreviations: n.s., not significant; YMRS, young mania rating scale. 
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Pathway analysis 

To detect whether a global expression pattern of a set of genes was 
significantly related to BD, pathway analysis was performed using 
Globaltest.^^ Gene set information was retrieved from the KEGG database 
(215 pathways; KEGG.db, version 2.5.0). The probe IDs were mapped to 
ENTREZ gene IDs using the Bioconductor annotation library Homo sapiens 
hugenelOstvl VI 2.0.0^^ for Molecular and Behavioral Neuroscience 
Institute Custom Chip Definition Files (MBNI CDF). 

Pathway scores were corrected for multiple testing using the Benjamini- 
-Hochberg method. To account for the size of the gene set, a probe label 
permutation test was performed.^ ^ The empirical P-value was calculated as 
the number of permutations obtaining a test statistics greater than, or 
equal to, the one observed in the original data, divided by the overall 
number of permutations applied. 

Detecting the associated genes in replicated pathways 
Any probe with a single probe P-value of 0.05 (as calculated by the 
Global Test) was considered to contribute to the significance of a pathway. 
The resulting list of probe IDs was mapped to ENTREZ gene IDs. To 
visualize potential crosstalk between significant pathways, a hierarchical 
clustering algorithm was used (average linkage method with the Jaccard 
coefficient as the distance measure) to group pathways with similar 
contributory gene content. An implicated SNP was defined as an SNP with 
a single SNP P-value of ^0.05. 

RESULTS 

Gene-expression differences in BD patients during euthymic and 
manic phases and in comparison with controls 
After correction for multiple testing, (i) 262 genes were 
differentially regulated in the comparison between euthymic 
patients and controls (177 downregulated (68%) and 85 
upregulated (32%) in BD euthymic); (ii) 216 genes were 
differentially regulated in the comparison between manic patients 
and controls (44 downregulated (20%) and 172 upregulated 
(80%)); and (iii) 22 genes were differentially regulated in the 
intrapatient comparison of euthymia and mania (all upregulated 
in BD manic; threshold: P < 2.59x 10~^). 

A set of 82 genes was dysregulated (40 upregulated, 42 
downregulated) in both BD euthymic versus controls and BD 
manic versus controls. A set of 16 genes was dysregulated in both 
the comparison of BD euthymic versus controls and BD euthymic 
versus BD manic (all upregulated; Figure 1). 

qRT-PCR 

qRT-PCR results for STABl show a trend towards the direction of 
mRNA changes observed in the microarray experiments 
(Supplementary Figure 1). In view of the small sample size, these 
do not surpass the significance level obtained through the 
performance of a f-test. 

Integration of gene-expression data with GWAS data 
As with genetic variation, genes that are differentially expressed in 
both the manic and the euthymic phases versus controls are not 
representative of a specific BD state but rather can be considered 
trait markers of the disease. Instead, these genes are more likely to 
be implicated in the etiology of BD per se. Therefore, the gene- 
based analyses focused on the 82 genes that were dysregulated in 
both the BD euthymic versus controls and the BD manic versus 
control comparisons. Data for 59 of these genes were available in 
the PGC-BD primary data set. Significantly, more genes were 
associated with BD than would be expected by chance (P = 0.04). 
A list of the top 10 associated genes, that is, those with a 
nominally significant gene-based P-value in the analysis of the full 
set of SNPs, is provided in Table 2. One gene — Stabilini 
{STABl) — had an empirical P-value of 1.9x10"^ which exceeds 
a Bonferroni-corrected threshold of P<6.1x10~'^ (that is, 
0.05/82). 

© 2014 Macmillan Publishers Limited 



Euthymia - Control Euthymia - Mania 




134 



Mania - Control 

Figure 1. Venn diagram. Number of differentially expressed genes 
across considered comparisons. 



Pathway analysis 

Analysis of KEGG pathways was performed for both the BD 
euthymic versus control and the BD manic versus control 
comparison groups. After correction for multiple testing, no 
pathway had a significant P-value in the BD euthymic versus 
control analysis. In contrast, the BD manic versus control analysis 
yielded nine significant pathways (Table 3). Clustering of 
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Table 2. List of the top 1 0 associated genes from GWAS of BD patients 
and controls 



Gene 


Chr 


nSNPs 


P 


minP 


STABl 


3 


71 


1.9x10""^ 


5.5x10"^ 


RBM14 


11 


23 


1.0x10"^ 


1.0x10""^ 


YWHAG 


7 


70 


1.2x10"^ 


8.3x10""^ 


SLPI 


20 


108 


1.3x10"^ 


8.2x10""^ 


HIST1H3H 


6 


60 


3.2x10"^ 


2.8x10"^ 


HNRNPC 


14 


107 


3.9x10"^ 


1.7x10""^ 


PCID2 


13 


90 


4.0x10"^ 


1.5x10"^ 


CDKN2D 


19 


47 


4.0x10"^ 


9.8x10""^ 


ANKRD36 


2 


18 


4.6x10"^ 


8.3x10"^ 


CHST14 


15 


52 


5.0x10"^ 


2.0x10"^ 



Abbreviations: BD, bipolar disorder; Chr, Chromosome; GWAS, genome- 
wide association study; minP, the original association P-value for the best 
SNP within the gene; P, the gene-based P-value considering the full set of 
SNPs; SNP, single-nucleotide polymorphisms. 



Table 3. Significantly dysregulated KEGG pathways in the comparison 
group mania vs controls 



Pathway ID KEGG pathway P No. of 

probes 



hsaOSOlO 


Alzheimer's disease 


2.5x10" 


3 


157 


hsa05012 


Parkinson's disease 


2.8x10" 


3 


113 


hsa05016 


Huntington's disease 


4.4x10" 


3 


171 


hsa04260 


Cardiac muscle contraction 


1.9x10" 


3 


70 


hsa00380 


Tryptophan metabolism 


1.9x10" 


3 


41 


hsa00190 


Oxidative phosphorylation 


2.6x10" 


3 


116 


hsa00140 


Steroid hormone biosynthesis 


3.9x10" 


4 


50 


hsa00980 


Metabolism of xenobiotics by 


8.7x10" 


4 


64 




cytochrome P450 








hsa03010 


Ribosome 


5.2x10" 


3 


75 
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" hsa00980:Metabolism of xenobiotics by cytochrome P450 



■ hsa00140:Steroid hormone biosynthesis 



hsa00380:Tryptophan metabolism 



hsa05016:Huntington's disease 



hsa05012:Parkinson's disease 



hsaOOl 90:Oxidative phosphorylation 



hsa05010:Alzheimer's disease 
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Height 

Figure 2. Mania associated KEGG pathways clustered on the basis of innplicated genes. CI, Cluster 1; C2, Cluster 2; C3, Cluster 3. 



significant pathways on the basis of gene-content sinnilarity 
generated three pathway clusters (Figure 2). Cluster 1 consisted of 
pathways classified as being involved in hunnan diseases, cluster 2 
consisted of nnetabolisnn pathways and cluster 3 contained the 
ribosonne pathway fronn the class of KEGG pathways with 
involvennent in genetic infornnation processing. Owing to the 
small nunnber of dysregulated genes, pathway analyses were not 
possible for the BD nnanic versus euthynnic connparison. 

DISCUSSION 

To identify new BD genes and pathways, the present study 
applied a novel, three-step convergent approach. First, candidate 
genes were identified through intrapatient connparisons of gene- 
expression profiles during euthynnic and nnanic states, and 
case-control connparisons. Second, differentially expressed genes 
were tested for association with the disease in existing BD GWAS 
data. Proof of a genetic association strengthens the hypothesis 
that a gene is innplicated in the etiology of the disorder, and that 
its differential expression is not merely a downstream epipheno- 
menon. Third, the implicated molecular mechanisms were 
investigated using pathway analytical methods. In total, 262 
genes were dysregulated in BD euthymic versus controls (85 
upregulated, 177 downregulated), and 216 genes were differen- 
tially expressed in BD manic versus controls (172 upregulated and 
44 downregulated. Gene-expression profiling identified 82 genes 
that were dysregulated in both the BD euthymic and BD manic 



versus control comparisons (40 upregulated, 42 downregulated). 
Intrapatient comparisons of the two mood states revealed 22 
dysregulated genes (all upregulated in euthymia). Of these, 16 
were downregulated in BD euthymic patients compared with 
controls (Figure 1). 

Whether euthymic BD patients can be viewed as 'healthy' is a 
matter of ongoing debate.^ In the present sample, euthymic 
patients differed from controls in terms of gene expression, and a 
large number of genes were dysregulated in patients irrespective 
of mood phase. Thus, our data underline differences between 
euthymic patients and controls, and suggest that similar, as well as 
specific, genes are dysregulated during euthymic and manic 
phases. 

In the second step, differentially regulated genes were tested 
for association in GWAS data using a gene-based approach. This 
analysis focused on genes that were dysregulated in both mood 
states versus controls. Of 10 nominally significant genes, STABl 
remained significant after correction for multiple testing (Table 2). 
Interestingly, STABl is located in close proximity to genes on 
chromosome 3p21.1 that have shown genome-wide significance 
in previous studies, that is, PBMR1 (-20 kb 3' of STABl) in a GWAS 
of mood disorder (P-3.6x 10~^),^^ and NEK4-ITIH1-ITIH3-ITIH4 
(~190kb 3' of STABl) in a meta-analysis combining BD and 
schizophrenia GWAS data (Praw = 8.4 x 1 0 ~ ^; Supplementary 
Figure 1).^^ To date, no specific gene has been conclusively 
shown to account for these association findings. The present 
results, however, suggest that STABl is a new and highly 
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Table 4. Genes in the pathways of interest 



Cluster 1 



Human disease 



Gene 



hsa04260 Cardiac muscle contraction 



hsaOSOlO Alzheimer's disease 



hsa00190 Oxidative phosphorylation 



hsa05012 



Cluster 2 



Parkinson's disease 



hsa05016 Huntington's disease 



Metabolism 



ATPIAh ATP1B3, ATP2A2, CACNA2D4, C0X4I1, C0X5A, C0X5B, C0X6B1, C0X6C, COX7A2, 
COX7A2L, C0X7B, C0X7C, C0X8A, CYCl, SLC9A6, TPM4, UQCRW, UQCRll, UQCRCl, 
UQCRC2 

ADAMW, ADAM! 7, APHIA, APP, ATP2A2, ATP5A 1, ATP5B, ATPSCl, ATP5D, ATP5E, ATP5F1, ATP5G2, 
ATP5G3, ATP5H, ATP5J, ATP50, BID, CALMI, CALM2, CALMS, CAPNl, CAPN2, CASP3, CDK5, CHP, 
C0X4I1, C0X5A, C0X5B, C0X6B1, C0X6C, COX7A2, COX7A2L, C0X7B, C0X7C, C0X8A, CYCl, FADD, 
FAS, GNAQ, HSD17BW, lUB, LPL, LRPl, MAPK3, NAEl, NCSTN, NDUFAl, NDUFA2, NDUFA3, 
NDUFA4, NDUFA5, NDUFA6, NDUFA7, NDUFA8, NDUFA9, NDUFBW, NDUFBll, NDUFB2, NDUFB3, 
NDUFB4, NDUFB5, NDUFB7, NDUFB8, NDUFCl, NDUFSl, NDUFS2, NDUFS3, NDUFS4, NDUFS5, 
NDUFS7, NDUFS8, PLCB3, PPP3CB, PPP3R1, PSENl, PSENEN, SDHB, TNF, TNFRSFIA, UQCRW, 
UQCRll, UQCRCl, UQCRC2 

ATP5A1, ATP5B, ATPSCl, ATP5D, ATP5E, ATP5F1, ATP5G2, ATP5G3, ATP5H, ATP5I, ATP5J, ATP5J2, 
ATP5L, ATP50, ATP6AP1, ATP6V0B, ATP6V0D1, ATP6V0E1, ATP6V1A, ATP6V1B2, ATP6V1C1, 
ATP6V1D, ATP6V1E1, ATP6V1G1, C0X17, C0X4I1, C0X5A, C0X5B, C0X6B1, C0X6C, COX7A2, 
COX7A2L, C0X7B, C0X7C, C0X8A, CYCl, NDUFAl, NDUFA2, NDUFA3, NDUFA4, NDUFA5, NDUFA6, 
NDUFA7, NDUFA8, NDUFA9, NDUFBIO, NDUFBll, NDUFB2, NDUFB3, NDUFB4, NDUFB5, NDUFB7, 
NDUFB8, NDUFCl, NDUFSl, NDUFS2, NDUFS3, NDUFS4, NDUFS5, NDUFS7, NDUFS8, PPAl, SDHB, 
TCIRGl, UQCRIO, UQCRll, UQCRCl, UQCRC2 

SEPTS, ATP5A1, ATP5B, ATPSCl, ATP5D, ATP5E, ATP5F1, ATP5G2, ATP5G3, ATP5H, ATP5J, 
ATP50, CASP3, C0X4I1, C0X5A, C0X5B, C0X6B1, C0X6C, COX7A2, COX7A2L, C0X7B, C0X7C, 
C0X8A, CYCl, HTRA2, NDUFAl, NDUFA2, NDUFA3, NDUFA4, NDUFA5, NDUFA6, NDUFA7, 
NDUFA8, NDUFA9, NDUFBIO, NDUFBll, NDUFB2, NDUFB3, NDUFB4, NDUFB5, NDUFB7, 
NDUFB8, NDUFCl, NDUFSl, NDUFS2, NDUFS3, NDUFS4, NDUFS5, NDUFS7, NDUFS8, PARK7, 
SDHB, SLC18A2, SLC25A5, UBA7, UBE2G1, UBE2J1, UBE2J2, UBE2L6, UQCRIO, UQCRll, UQCRCl, 
UQCRC2, VDAC3 

AP2S1, ATP5A1, ATP5B, ATPSCl, ATP5D, ATP5E, ATP5F1, ATP5G2, ATP5G3, ATP5H, ATP5J, 
ATP50, CASP3, CLTB, CLTC, C0X4I1, C0X5A, C0X5B, C0X6B1, C0X6C, COX7A2, COX7A2L, 
C0X7B, C0X7C, C0X8A, CYCl, DCTNl, DCTN2, DCTN4, DLG4, DNAL4, GNAQ, HDAC2, NDUFAl, 
NDUFA2, NDUFA3, NDUFA4, NDUFA5, NDUFA6, NDUFA7, NDUFA8, NDUFA9, NDUFBIO, 
NDUFBll, NDUFB2, NDUFB3, NDUFB4, NDUFB5, NDUFB7, NDUFB8, NDUFCl, NDUFSl, 
NDUFS2, NDUFS3, NDUFS4, NDUFS5, NDUFS7, NDUFS8, PLCB3, P0LR2B, P0LR2E, P0LR2H, 
P0LR2K, P0LR2L, PPARG, REST, SDHB, SLC25A5, SODl, UCPl, UQCRIO, UQCRll, UQCRCl, 
UQCRC2, VDAC3 

Gene 



hsaOOMO Steroid hormone biosynthesis 
hsa00980 Metabolism of xenobiotics by 

cytochrome P450 
hsa00380 Tryptophan metabolism 

Cluster 3 Genetic information processing 



COMT, CYPIBI, HSD17B12, HSD17B7, SRD5A1, SRD5A3, UGT2A1, UGT2B15 

ADH5, LDH3B1, CYPIBI, CYP2E1, GSTKl, GSTM5, GSTPl, MGSTl, MGST2, MGST3, UGT2A1, 

UGT2B15 

ACATl, ACAT2, AFMID, CYPIBI, ECHSl, GCDH, HADHA, MAOA, MAOB, OGDH, TPHl, WARS 
Gene 



hsa03010 



Ribosome 



FAU, RPLIO, RPLIOA, RPLll, RPL12, RPL15, RPL17, RPL18, RPL19, RPL23, RPL24, RPL26, RPL27, 
RPL28, RPL29, RPL3, RPL30, RPL31, RPL34, RPL35, RPL35A, RPL36AL, RPL37A, RPL38, RPL4, RPL41, 
RPL5, RPL7A, RPL8, RPLPO, RPLPl, RPLP2, RPSl 1, RPS13, RPS15, RPS15A, RPS16, RPS20, RPS23, 
RPS24, RPS25, RPS27L, RPS4X, RPS4Y1, RPS5, RPS6, RPS7, RPS9, RSL24D1, UBA52 



pronnlslng candidate gene in this region. A recent pathway 
analysis of a schizophrenia GWAS, which identified STABl as one 
of 15 candidate genes, provides further support for this hypo- 
thesis.^° Possible approaches to the further investigation of this 
hypothesis are genetic fine nnapping of the associated regions, 
and studies of effects of the associated SNPs on STABl expression. 

STABl encodes a nnultifunctional scavenger receptor, and its 
expression is induced during chronic inflannnnation and tunnor 
progression.^^ Post-nnortenn gene-expression studies in the 
dorsolateral prefrontal cortex have revealed differential expression 
of STABl in BD and schizophrenia patients connpared with controls 
(see www.stanleygenonnics.org), although the direction of these 
changes were inconsistent across studies. Enignna, which is the 
largest consortiunn conducting nneta-analyses of neuroinnaging 
and genetic data worldwide, perfornned nneta-analyses of 
structural neuroinnaging data, and reported one significant 
association in the above nnentioned 3p21.1 region, that is, 
between average caudate volunne and an SNP in STABl 
(P = 7x10~^ Supplennentary Figure 2).^^ The caudate nucleus is 
part of the ventral striatunn, a brain region innplicated in BD.^^ 



STABl binds and internalizes SPARC (secreted protein acidic and 
rich in cysteine), a protein which is also involved in 
inflannnnation.^'^ SPARC has been innplicated in the regulation of 
depression and affect anxiety-related behavior in nnice.^^ As STABl 
was upregulated in both the BD euthynnic- and the BD nnanic 
versus control connparisons (fold change =1.5; ^=1.3x10""^ and 
P = 2.5x10-4, respectively), this gene could be a nnarker for BD 
per se, irrespective of nnood phase. Further studies are warranted 
to elucidate the role of STABl and its interplay with SPARC in BD. 
Mouse nnodels with STABl deficiency are available^^ and represent 
a potential tool for further analyses. 

Pathway analyses of genes fronn the BD nnanic versus control 
connparison revealed nine pathways with statistically significant 
differential expression (Table 3). Clustering of these pathways on 
the basis of gene-content sinnilarity generated three clusters 
(Figure 2). These clusters were nnainly driven by dysregulated 
genes whose biological functions have been innplicated in BD by 
previous research. 

In cluster 1, enrichnnent of genes involved in energy nnetabo- 
lisnn and the nnitochondrial connplex (Table 4) was observed. 
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Transcriptional studies of multiple brain regions have implicated 
these genes in BD (for review see Konradi et al}^). Mitochondrial 
dysfunction is an important component of the pathophysiology of 
BD,^^ and mitochondrial modulators are considered potential 
therapeutic targets.^^ The cluster 1 pathways were also enriched 
in ubiquitin proteasome genes (Table 4). Dysregulation of these 
genes has been reported in brain- and blood-based gene- 
expression studies of patients with psychosis.^^ Cluster 1 pathways 
also contained a number of genes involved in inflammation 
(Table 4). Cytokines have been implicated in BD and other 
psychiatric disorders, although these data are limited (for example, 
Guloksuz et al.^^). Finally, genes from the pathways in cluster 1 
have been implicated in protein trafficking (for example, clathrin), 
and protein trafficking abnormalities have been linked to the 
pathophysiology of BD (for a review, Schubert et al.^^). 

In the cluster 2 pathways, differential expression was mainly 
driven by cytochrome P450 proteins, which catalyze many of the 
reactions involved in drug metabolism and the synthesis of 
lipids.^^ Cluster 2 also contained proteins involved in the 
biosynthesis or degradation of neurotransmitters that have been 
implicated in BD previously (Table 4).^^'^"^ 

The cluster 3 pathway consisted of ribosomal genes. The 
function of the ribosome is protein synthesis, and research 
suggests that disturbances in protein synthesis may be related to 
psychiatric disorders. For example, the control of protein synthesis 
via elF2-alpha kinases is activated by memory processes, as 
mediated by glutamate receptors, and many of the genes in this 
network have been associated with BD.^^'^^ In this context, it is of 
note that during the mania phase, 80% of the differentially 
regulated genes were upregulated, whereas during the euthymia 
phase, 68% were downregulated in the present patients. BD 
patients may become more physically active during a manic 
phase, and comparable changes in gene expression have been 
found in studies of exercise and training.^^'^^ Even if dysregulation 
of these 'exercise genes' is only a marker for dysregulated 
processes during mania and does not represent a causal 
mechanism, it is nonetheless of potential value in a clinical 
context, as early awareness of such dysregulation, together with 
prompt and effective counter strategies, could prevent the 
development of a fully-fledged manic episode.^^'^° Combining 
the assessment of clinical symptoms and biological biomarkers 
could enhance the specificity and selectivity of mood-change 
prediction. 

Although differential pathway expression was found in the BD 
manic versus control group, no pathway showed association in 
the BD euthymic versus control group. One reason for this may be 
that euthymia is a more heterogeneous state than mania. Thus, 
during mania, when patients exhibit very pronounced and severe 
symptoms, specific dysregulated genes are pronounced. As manic 
and euthymic gene expressions were compared on an intrapatient 
basis, the results are not attributable to inter-individual 
differences. 

A major limitation of the present study was a lack of knowledge 
concerning the influence of medication-related factors. Thus, the 
possibility that the differences observed between patients and 
controls were because medication cannot be excluded. In 
particular, some of the alterations discussed above are influenced 
by lithium,"^^"^^ which was taken by 50% of the present patients in 
both the manic and the euthymic phases. As the direction of gene 
dysregulation substantially differs between the two states, it is 
unlikely that these differences were entirely caused by medication. 
Although no significant differences were observed in terms of the 
classes of medication prescribed during the manic and euthymic 
phases, the dosages differed significantly between the two mood 
phases. Therefore, the possibility that differential gene expression 
was partly because differences in medication dosages cannot be 
excluded. Nevertheless, the findings of the parallel GWAS 

Translational Psychiatry (2014), 1 -7 



approach, which is trait-specific and robust against such medica- 
tion effects, support our results. 

Further limitations were the small sample size and the fact that 
the patients were not stringently assessed with respect to time 
points, medications, and alcohol and tobacco use. In addition, 
patients did not have the same scores in their manic and euthymic 
phases but were selected in a top down manner with respect to 
the highest difference in YMRS scores. All of these factors 
increased the heterogeneity of the sample and — given the 
stringent threshold for significance — decreased the chances of 
identifying significant gene-expression differences. 

Another limitation is that gene expression was analyzed in 
peripheral tissue, which does not entirely reflect gene-expression 
changes in the brain.^^ However, many of the genes that were 
differentially expressed are also found in the brain.^^ Furthermore, 
previous authors have demonstrated the value of analyzing 
peripheral tissue to explore molecular mechanisms and detect 
new genes and pathways for BD.^'"^^ Finally, intrapatient compar- 
isons of mania and euthymia were only possible using peripheral 
tissue, and all potential biomarkers derived from transcriptional 
profiling must be measurable in easily accessible tissue if they are 
to be developed for routine clinical use. 

In summary, the present study applied a novel three-step 
convergent approach that involved intrapatient and case-control 
comparisons of euthymic and manic gene-expression profiles. To 
our knowledge, this is the first study to compare gene-expression 
profiles and dysregulated pathways in this manner. The findings 
suggest that besides differing from controls during their manic 
phases, BD patients differ from controls during euthymic phases. 
The results provide support for previously implicated molecular 
mechanisms, and suggest new BD genes and mechanisms. We 
hypothesize that genes and pathways implicated in the manic 
phase relate to the mania state rather than to disease-specific BD 
traits. Differentiation between trait- and state-specific gene 
expression through intrapatient comparisons may be of value in 
identifying novel vulnerability genes and medication targets, and 
the delineation of biomarkers for the early detection of mania 
onset. Our data suggest STABl as a new and highly promising 
candidate gene for BD. As STABl was upregulated during mania 
and euthymia, further investigation is warranted to determine 
whether this gene is a potential medication target or biomarker 
for BD. 
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